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Abstract 

The effective action is derived for a self-interacting theory with a 
finite fixed 0(2) charge at finite temperature in curved spacetime. We 
obtain the high temperature expansion of the effective action in the 
weak coupling limit. In the relativistic temperature, we discuss about 
the phase transition in a homogeneous spacetime. 
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I. INTRODUCTION 



Since the discovery of the Bose-Einstein condensation (BEC) much ef- 
fort has been made to understand the effect of the non-vanishing chemical 
potential on thermodynamical quantities in quantum field theories in curved 
spacetime, where the spatial section may have a boundary. The high tem- 
perature behaviors of the grand thermodynamic potential in an arbitrary 
static spacetime have been investigated by utilizing the heat kernel and 
zeta-function regularization technique [Q]. Especially, for an Einstein uni- 
verse which is a homogeneous space and have no boundary, the exact results 
for the thermodynamical quantities are derived in the non-relativistic limit 
and in the relativistic limit 0. Recently D. J. Toms studied the BEC in a 
static spacetime with a possible spatial boundary and derived the universal 
character-the portion of the charge in the ground state to the total charge 
is identical to the result in the flat spacetime ||. But to our knowledge all 
of them does not treat the interaction term. 

For the system having the interaction terms, the effective actions ( poten- 
tial ) have been only calculated in the flat spacetime. Kapusta, by studying 
the effective equation of motion for fields, showed that charge conservation 
can affect the phase transition of a theory [Q. Haber and Weldon studied the 
effects of a net background charge on ideal and interacting relativistic Bose 
gases in the large iV approximation ||. More recently Benson, Bernstien, 
and Dodelson exactly calculated the one-loop finite temperature effective 
potential for a self-interacting theory with a fixed 0(2) charge 0. 

In this paper we calculate the effective action for a model with an in- 
teraction term and a fixed charge at finite temperature in curved spacetime, 
and study the effects of the interaction on the phase transition. The simple 
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model under investigation, is described by the matter Lagrangian density 

- £ = -(^i^i + Vfadpfa + m 2 2 ) + W 2 + ^\ (1) 

where £ is a numerical factor and R is the Ricci scalar curvature. This model 
is 0(2) invariant, and there is a conserved charge Q associated with 0(2) 
global symmetry. 

In Sec. II we set up the formalism for evaluating the effective action at 
finite temperature and charge. To evaluate the effective action, we introduce 
the Riemann nomal coordinates and use the heat-kernel. In Sec. Ill by 
using the dimensional regularization method, we derive the high temperature 
expansion of the effective action in the case that the coupling constant A is 
small, and briefly discuss about the renormalization. In Sec. IV using the 
result of Sec. Ill, we discuss about the phase transition of the theory in the 
relativistic temperature limit in a homogeneous spacetime with no boundary. 
Finally we summary the results. 

The spacetime we discuss is the ultrastatic spacetime M. = R x E. 

ds 2 = dr 2 + g ij (x)dx i dx i , (2) 

where the Wick rotation r = ix° has to be understood. 

II. ONE-LOOP EFFECTIVE ACTION 

In this section we compute the effective action at finite temperature and 
finite charge. According to Ref.||, §, the grand partition function Z[{3, /i] = 
e -/3fi(/3,/i)^ Q({]^ij) being the thermodynamic potential, can be expressed as a 
path integral: 

Z = N( [d^d^expi-S^i]), (3) 
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where N is a constant and the integration has to be taken over all fields 
0(r, x) with periodicity ft = k with respect to r. The action S is 



S 



dr J d x^fg 
X 



{Vfaf + (V0 2 ) 2 



2 J,2 



m (j) 



V 2 ^ 



HRfa] + || fa + ¥{fafa ~ Mi) - ^ 



(4) 



where <jr = <fi{ + (f)^, fa = and fi is the Lagrangian multipler related to 
the conserved charge Q. It is the same form as the flat spacetime. 

In order to calculate the effective action, some perturbative approach 
must be adopted. The usual one is the loop expansion J7|. We now expand 
the action around a background field configuration, fa(x) = fa(x). fa are 
determined by the equation J£ ^. = ^. = , where T is the effective action. In 
general fa(x) is not a constant field in curved spacetime. Then the action in 
Eq.([|) is expanded in powers of fa(x) = fa(x) — fa(x): 



s = s<°> + + s< 2 > + 



(5) 



The zeroth order term is just the action evaluated at fa(x): 
= p f d 3 x^ 



\[{Vfaf + (V0 2 ) 2 + m 2 2 + m 2 } + ^ ~ ^ 2 



where 



(6) 

L >\ + fa\. The first order in Eq.(Kp, which is linear in fa(x), can 



be neglected. The second order term in the action is 

S {2] = i / d*x^g [ d*y^g oXrl.M,, </(//). 



(7) 



where the 2x2 matrix A4ij is the second functional derivative of the action 



with respect to the field evaluated at fa, 

62 S 

Mi 



Sfa(x)Sfa(y) 



(8) 
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Since the boundary condition at finite temperature is periodic for a bosonic 
field, we may expand <^(t, x) in Fourier modes: 

1 oo 

^{r,x) = -j= ( n We M " T (9) 



-oc 



with uj n = ^p. Substituting Eq.(g) into Eq.([7|), we obtain 



Z {1) = [#'J[4, J exp 



(10) 



where .M \j = (— A 3 + + a> 2 — /i 2 — 2fiuj n eij). Here A 3 is the Laplacian on 
the spatial section S. V^- and are given by 

Vij = £R8ij + fhl(4>)4>i^j + rh 2 2 (4>)(8ij - <fo^), e {j = \\ q ) ' 
where 

m 2 M = m 2 + ^ 2 , (12) 
m 2 M = m 2 + ^0 2 , (13) 

(14) 



In terms of a loop expansion the effective action is then given by 

r[0] = ^°)[0] + r one _ ioop + ---, (is) 

where the one-loop effective action T on e-ioo P is given by 

F one-loop = ~ In Z {1) . (16) 

It is note that the functional at finite temperature in d = 4 dimensions 
is completely reduced to 3— dimensional expression. Using this fact Eq.(|T6|) 
can be calculated by the proper time formalism in 3-dimensions. 



5 



Following DeWitt ||, let us write 

i oo f poo A c 

lnZ« = 5 E / -Me" 
= ^Jj 3 x^gT 



sH„ 



x) 



with 



(17) 



(18) 



-A 3 + V t j +L0 n - fX - 2[XU) n €ij. 

Since the eigenvalues of the operator H n are not exactly known for a general 
background metric g, it is necessary to make use of some approximation 



scheme. The useful approximations are the Riemann normal coordinates [10 



and the weak field expansion. We will use the result of Ref.yjJ. Following 
the ansatz for the heat kernel suggested by Jack and Parker, we use the 
nonlocal form e~ a ^ v ~^ R ^ in the heat kernel |I2] . The explicit form of T, up 
to adiabatic order four, is given by ( see Appendix A ) 

d N p 



^1 g roods 

ft n=oo •'O ^ 



(2n) 



N 



% + h + I 2 + h + h + -- ■) 



r + r 1 + r 2 + r 3 + r 4 + 



where 

h 
h 
h 

h 

h 



e -s{ P 2 +V-\R) 

o, 

_ se -s{p*+V-lR) 
_ se -s( p 2 + V-lR) 



\R-2A aP {8' 
o 



a/3 



sp a p 13 ) 



sV aP a - 3v4 3 a/3(-^ a V + -S 2 pW) 



(19) 

(20) 
(21) 
(22) 

(23) 



1 ~ 
2 



_|_ ^7^/35 + £ C 
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li?-2A 2a/3 (V-^pV) 



: R-2A 2a(3 S afS 



6 
1 

6' 



+ 2sA 2a /3 
2 
9 ; 



(24) 



Here V = Vij + u\ — y 2 — 2fMj n eij. After integrating about the momentum 
p, we get 

r°° ds 1 



r ° = tr / -.Ett^ 1 
ri = r 



/5 n (47TS)- 

r 3 = o, 



(47TS) 



— ny 

6 



30 180 ^ 7 180 ^ 



(25) 
(26) 

(27) 



where R, R a p-yS and R a p are curvatures on the 3-dimensional hypersurface 
£ The evaluation of the T and T 4 is not easy because the potential V is 
not a diagonal matrix. However if we use the property of the trace, we can 
diagonalized e~ s ^ v ~^ (see Appendix B). If we choose the matrix S as 



S 



(ml - m 2 2 )<\)x<\>2 - 2fiu n (ml — m\)4>i4> 2 — 2fiuj r , 
X+-C A_ - C 



where 



C=(£--)R + m\4>{ + mi(/>i + u, 



(28) 



(29) 



and A± are the eigenvalues of the matrix V — 



A ± = (e-^ + m 2 + ^ 2 + ^-^±i 



- 16A£ 



(30) 



7 



Then we can rewritten the Tq and as a more tractable form: 



1 



r = j}U*«r m , - s ' N/2 ^ 

P n J C 



r 4 = ^u^r N/2 

P n 



ds 
s 

00 ds 



-N/2+2 tl 







e -sA+ 


















e -sX 



(31) 



— S~ 1 avs + d 2 
6 



/3 n JO S [ 



-sAj 



+ e- x -)Ca 2 -±n^) 



18 



A_ 
18 



A 



+ -A0102 D (0102) 



-sAj 



-sA_^ 



(32) 



where d 2 = ±DR + ^R^R^ - -L.R aj3 R a P. 

Now, let us consider limiting case. First of all, we consider the case that 
the metric is flat. In flat spacetime the background fields is a constant 
because the spacetime is homogeneous. Then only T survive. Using the 
formula 



G ri 



dT 

dm 2 



-^/^/;^- + e-^-)) 



(2n) N 

d N p ( 1 1_ 

(2^p ly + A+ + p 2 + A 



(2n) N dm 2 
l r A 2 



In 



w n -A 2 -^- 1 ^ 2 ^ 



(33) 



we obtain, after sum over n, 



(3 J (2tt) 



AT 



+ 2 ln(l + e-^+) 
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+ 2 ln(l + e 



(34) 



where E 2 ± = p 2 + M] + y 2 ± [V(p 2 + M 2 ) 1 / 2 + f^ 4 ] 1/2 with M 2 = m 2 (l + 
Y^cj) 2 ). In this case, the effective action ( potential) is given by 



+ ln(l + e~ pE+ ) + ln(l + e~ pE - 



(3{E + + 



(35) 



This is exactly the same form as the previous result in the flat spacetime M. 

Second, for the case A — > 0, \± — > (£ — |)i? + m 2 + (cj n ± ifi) 2 . Then 
T and T 4 go to the previous result because OV — > 0. For /i = 0, A+ = 



(f - §)# + m 2 + io 2 n + ±A0 2 and A_ = (f - 



u 2 + \\<p 2 . In this 



case the effective action goes to the previously known result [|TJ. 

III. HIGH-TEMPERATURE EXPANSION 



In this section we will calculate T and in the high temperature limit 
and the coupling constant A is small. The regularization method is the 



dimensional one |13[ . After the regularization, we take N — > 3. Fisrt of all, 
let us calculate IV For small A, Tq can be written as 



1 00 N 

\ E (4-)-^ 2 r(-f ; 

P n=-oo Z 



+ 



M 2 + co 2 n -^ 2 - 



,(n=0) 



+ r« + rr + ■ 



(2) 



- 16^ 2 w 2 



1 / A 2 - 



1 - 16/^ 2 



(36) 



where 



r 



(n=0) 



M 



Ai 2 + 



A TO 
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r (l) 

1 



F (2) 
1 



+ 



m 2 - ^ - h 2 

o 



(37) 



1 AT CO s 

i(4 7 r)^/ 2 r(-|)2E{K + 



/3 
+ 
1 



M 2 + (u n + ifi)' 



2 2 



(38) 



/3 (4 7 r)-^ 2 r(-^)2f 



.iV A 2 4 

i 

' 2 144/xu; ri 



n=l 



M 2 + (cu n - i/x) ; 



; iV A 2 4 
' 2 144/za; n 



M 2 + (w n + i/x) 1 



(39) 



where M 2 = m 2 — /s 2 + £R+ ^(j) 2 . The evaluations of the Matsubara sums in 
Eqs.( p7| , p8| , p9D are not difficult |I4|. By using the dimensional regularization, 
we obtain 



,(0) 



,(2) 



1 1 

/?6tt 
V 1 



(M 2 - /i 2 + -0 2 ) 3 / 2 + (M 2 - /i 2 - -0 2 ) 3 / 2 
6 6 



(40) 



45 f3 4 12/3 2 
2 1 
3 



(M 2 - 2^ 2 

2 



1 1 , J 2 M 2 . 

- + - 7 + In F- , 

e 2 47r 



M 4 
(4^ 



7T 



3 
2' 



(^ n M^) 2fc C(-3 + 2fc + 2r) 
A 2 4 



(41) 



288tt 2 



11 J 2 M 2 . 
e 2 47r 



1 



+ £ C(--,Mr + l) 

fc,r>0,fc+r^0 



( _ ir( ^ ) 2, ( ^ ) 2, c(2fc + 2r+1) 

where e = 4 — d and 7 is the Euler constant 0.577..., and 



(42) 



_ N _ k+r T(s-% + k)T(2s-N + 2k + 2r) 
1 2' ' ] 1 j A;!r!r( S -f)r(2s-iV + 2A;) 



(43) 
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Now let us calculate the T 4 . After some algebra, we get 



+ 



87T/3 



471-2 



1 1 . ,/3 2 M 2 , 

7 + 2 7 + ln( ^ 



+ J2 C{±k,2r) 

k,r>0,k+r^0 



,ijj£ 2 (3M 2k 
■ 2vr ' 1 2vr 



A 



C(l + 2k + 2r) (o2 - — n<f> 



M 2 — /i 2 + ^ 2 
o 



+ 



M 2 - f - -J 2 
o 



A 7- 

(oa - — □</>' 



M 2 — fj 2 + -J 2 
o 



m 2 -v - h 2 

o 



A 2 



288?r 2 



6 m 



)□( 



E c(i Mr + i)(-i)*(£r 



fc,r>0 



^) + 30i02n(0i0 2 ) 
^)*(^C(3 + 21fc + 2r), (44) 



where we have calculated up to A 2 order except n = term. The final result 
for the effective action is given by 



r[0] = s[4>] - -p jf d 3 x^ (if= 0) + r« + rf + r 4 + • • ■) (45) 

up to A 2 order. 

Let us consider the infinite part of the effective action. It is given by 

M 2 A 



^infinUe[4>] = ~ d 3 Xy/g2 



1 



(47T) 



"7 + r + In W' 



- + ^ 4 + («2--n0 2 ) 



(46) 



where /i is a arbitrary mass scale. It is note that the value of the pole part 
differs from one that calculated by the zeta function regularization method [|T], 



TM. In the case A — > 0, this reduce to the known result. These infinite terms 
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can be removed by absorbing it into the gravitational Lagrangian and by 
introducing the counter terms | 15| . Up to terms which are total divergence, 



the coupling constant and mass counter terms are only needed. This is 
coincide with the fact that at one loop level, there is no wave functional 



counter term [H3] . 



IV. BOSE-EINSTEIN CONDENSATION 

In quantum field theory, BEC is interpreted as a symmetry breaking 
effect in flat as well as in curved spacetime || ^ ||. In order to discuss 
the BEC in the case of a relativistic Bose gas, it is important to clarify in 
what sense we mean that the temperature is relativistic. In flat spacetime, 
it means that T ^> M. In curved spacetime, in addition to T ^> M , we 
require T> l-R] 1 ^ 2 , where \R\ is the magnitude of a typical curvature of the 
spacetime ||. 

The high temperature effective action, from Eq. fl4"5|) , has the form in the 
relativistic temperature, 



T[4>] = S[4>] + T 



one— loop 



7T 2 1 1 9 9n 

i + T^( M 2 -V) + 



45 p A 12P 2 



(47) 



At first we discuss about the case that the spacetime is homogeneous. The 
translational invariance implies that the background field is a constant ( This 
is always not true, For A = 0, see Ref. 0, [IT])- 111 this case, the equation of 
motion for fa implies 

(m 2 - /i 2 + iR + ^ 2 + ±T 2 )fa = 0. (48) 
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Thus the effective action have two minima: 



6 
A 
6 
A 



V 2 -(m 2 + ZR+-T 2 



/i 2 -m 2 (T) 



0, with unbroken symmetry, 
A 



(49) 



with broken symmetry. This is the same form as the result in the flat space- 
time except for changing m 2 — > m 2 + £R §. The critical temperature when 
phase transition occurs is 

T 2 = ^[n 2 (T)-(m 2 + ZR)} (50) 

In general, // 2 has the temperature dependence due to the charge conserva- 
tion. 

Now consider the expectation value of the charge operator Q which is 
given in terms of the effective action by 



Q 



1 OT 



(51) 



in the high temperature approximation. Qo is the charge in the ground state 



and Qth is the charge in the thermal excited state. If T is high, 



0. Then 



Q = ^T 2 V, where V is the volume of the spatial section E. As T decreases, fi 



must increase until the temperature becomes T c . From Eq.( 
the critical temperature is rewritten as 

/3Qx x ' 2 



and Eq.(|5 



\ V 



{m 2 {T))- l f 2 . 



For T < T c , it is easily seen that 

Qo = Q 



l - 



(52) 



(53) 
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which is identical form to the result in flat spacetime ||. 

Up to date we have restricted the value of the fields (pi to be a constant. 
But in curved spacetime or in the case that there are boundaries, as others 
showed, the value of the fields (pi is not a costant even for A = 0. In addition, 
the critical temperature differs from our result because there is a non-zero 
lowest eigenvalue for the Laplacian operator 0, O . 

In an inhomogeneous spacetime, the equation of motion for the fields (pi 
is given by 

□0. _ ( m 2 _^ + iR+ 1 02 + A T % = 0. (54) 

This is a nonlinear equation, so it is not easy to find the analytic solution. 
Therefore it is difficult to discuss about the phase transition of the system. 

V. DISCUSSION 

In this paper we have calculated the effective action for the 0(2) invariant 
model with quartic self-interaction and global 0(2) fixed charge at finite tem- 
perature in curved spacetime. We calculated the effective action by using the 
modified heat kernel technique suggested by Jack and Parker and introducing 
the normal coordinates expansion method. We also used the fact that a 2 x 2 
matrix can be diagonalized in general. We showed that the coefficient a± - 
the Minakshiuundaram-Seeley- DeWitt coefficient- does not appear similarly 
to the result at zero temperature case. 

Using the effective action, we discussed about the phase transition in the 
homogeneous spacetime. It is shown that the charge condensation property 
is the same as the flat spacetime case, except for the changing of the mass. 
The form of the condensated charge fraction to the total charge is the same 
as the non- interacting case (Eq. (|53|) . For the inhomogeneous spacetime we 
have not discussed about the phase transition because the equation of motion 
is nonlinear. 
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Even if we have studied the simple 0(N = 2) model, our result can be 
generalized to an arbitrary potential type and to the N > 2 model. But for 
N > 2 one will be confronted with the diagonalization problem of the matrix. 

For conformally ultrastatic spacetime the effective action may also be 
derived by using the conformal transformation technique. 
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APPENDIX A 



In normal coordinates, an operator 0(x, x') may be expanded in powers of y 



about the point x', where y 



fj, = dx»(r) 
dr 



T=0 



, x{j) being a geodesic connecting 



x and x' [TD[. Then, up to adiabatic order four 



where 



c 2 

Ci 



-a + V) = C 1 + C 2 + C 3 + C 4 + - 



(-d 2 + V-~R), 
6 

-(A 2a ^ + B^d^d v - ±J2), 
-(A 3a p iy a yPyW 2 + B^ aP7 d„y a yPyW v - V a y a ) 
-(A 2a/3 B^ s y a y^d^y 5 d u + A w y a y^y s d 2 



1 ~ 



+ B^ s d,y a y^y>d v - -V^y a y^ 



(55) 

(56) 

(57) 
(58) 

(59) 



with 



.4 



2a/3 



:R, 



a/3i 



(60) 
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^4q/3 7 5 = -^Ra/3-.-yS + ^^/^V + Ij^^aP^fS, (62) 

^ = ~ (63) 

b^,, = \R":^'Y2 6 " VRa ^ (64) 

1 „ 1 „„ 1 



+ ( 2Q^V/3: 7 5 + Jc-R^aXpR^S ~ ^ R ^ R Vs) ( 65 ) 

By using the Schwinger expansions 

e -(Ho+ff I ) s = e -H 0S + (_ s ) f 1 due'^-^'Hte-^ " (66) 

Jo 

+ (-s) 2 / c/wu / dve-^^Hje-^-^'H^-^ 1103 + ■■■, 
Jo Jo 

one may evaluate the heat kernel 

K = e (- D +^ 

Jo 

■(C 2 + C 3 + C 4 )e-" ( - s2+ ^^ )s 

+ (_ s )2 Z -1 ^ / 1 dve -(i-u)(-a a +^-|«)- C r je -«(i-«)(-a a +v-|H)- 
jo Jo 

■C 2 e-™ ( -° 2+ ^-^ )s + ■ ■ • . (67) 
In local momentum space, for any operator 0(x, x') 
(x\0(x,x')\x') = (y\O(y,x')\0) 

= jdpdq(y\p)(p\0(i^-ix')\q){q\x), (68) 

where y — x — x'. Then, using the trace property tr(ABC) = tr(CAB), one 
can obtain 

r d N v 

ti(x\K\x) = try j^£j(Io + Ii + h + h + h + .-.)• (69) 
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APPENDIX B 



For any real 2x2 matrix A, a matrix S can be chosen so that S 1 AS has 
the diagonal form. Let 



a b 
c d 



A 



a 



v7, 



a 



v7y 



A 



A_ 



ft 



(70) 



where X± are eigenvalues, and (^j and (^j are the corresponding eigenvectors. 
They are given by 



2 

::) ■ ' 



[a + d)± ((a + d) 2 - 4(ad - be)) 



1/2 



' b N 
, A_l — a , 



-5, 



' b " 
,A_ - a ; 



If we choose the matrix S 1 as 



5 = 



then 



S^AS 

From this we can show that 
S~ l e tA S 



a (3 \ 
7 5 J' 

/A + 
A_ 



e 5 lA5 = exp 



X + t \ 
A_ ) " 

tre A * = tr(5- 1 e^5) = (e A +* + e A -*), 

tre At B = tr{S- 1 e At SS- 1 BS) 
e x+t 



3 A +* 



tr 







(71) 
(72) 



(73) 
(74) 

(75) 
(76) 

(77) 
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